Dwell time symmetry in random walks and molecular motors 
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The statistics of steps and dwell times in reversible molecular motors differ from those of cycle 
completion in enzyme kinetics. The reason is that a step is only one of several transitions in the 
mechanochemical cycle. As a result, theoretical results for cycle completion in enzyme kinetics do 
not apply to stepping data. To allow correct parameter estimation, and to guide data analysis and 
experiment design, a theoretical treatment is needed that takes this observation into account. In this 
paper, we model the distribution of dwell times and number of forward and backward steps using 
first passage processes, based on the assumption that forward and backward steps correspond to 
different directions of the same transition. We extend recent results for systems with a single cycle 
and consider the full dwell time distributions as well as models with multiple pathways, detectable 
substeps, and detachments. Our main results are a symmetry relation for the dwell time distributions 
in reversible motors, and a relation between certain relative step frequencies and the free energy per 
cycle. We demonstrate our results by analyzing recent stepping data for a bacterial flagellar motor, 
and discuss the implications for the efficiency and reversibility of the force-generating subunits. 
Key words: motor proteins; single molecule kinetics; enzyme kinetics; flagellar motor; Markov 
process; non-equilibrium fluctuations. 
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Introduction 

Progress in single molecule techniques has enabled ob- 
servations of single steps in many motor proteins, and 
accurate measurement of the distribution of dwell times, 
i.e., the periods of little or no motion between steps. Ex- 
amples are the forward and backward steps of proces- 
sive molecular motors like kinesin [3,0, S 13, myosin 
V [1 H 0, H, cytoplasmic dynein | or RNA poly- 
merase II Oil : and stepwise rotations in ATP synthase 

[O, [H, M, [II, [H, El, [13 and the flagellar motor [lj]. 
Observations of steps and dwell time distributions offer a 
route to gain insight into the microscopic mechanism and 
detailed motion of such systems, beyond what is available 
through knowledge of average turnover rates alone. 

Dwell times are examples of first passage times, which 
have been extensively studied in the theory of random 
walks [l9| . An important modeling step is therefore to 
formulate a first passage problem that describes the ex- 
perimental situation. Close examination of step trajecto- 
ries with high time resolution from several motor proteins 

fllllilHJJ^EEQElEEIll re- 
veals that the steps are very rapid events compared to 
typical dwell times. In a discrete-state description, it is 
therefore reasonable to identify a step with a single tran- 
sition, and the direction of the step with the direction of 
that transition [2(| Hl| . The identification of steps with 
single transitions follows naturally from the assumption 
that each state has a well-defined average position, and 
is a basic assumption in this paper. 

Motor proteins are cyclic enzymes, but contrary to ear- 
lier assumptions [IH, HI] , the steps and waiting times ob- 
tained in many single-molecule experiments arc not well 
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described in terms of the cycle completions in enzyme ki- 
netics. This was recently demonstrated for models where 
a single cycle accounts for both forward and backward 
steps (2Cj. In this case, the average number of forward 
and backward steps differ from the average number of 
completed forward and backward cycles, and the average 
dwell times between steps differ from the average cycle 
completion times. The differences can be significant even 
in conditions where almost no backward steps occur. The 
basic reason for this difference is that a motor with one or 
more intermediate states per cycle is in a different state 
immediately after forward and backward steps. There- 
fore, an experimental trajectory of forward and backward 
steps does not give explicit information about completed 
forward and backward cycles. 

The different states obtained just after forward and 
backward steps also implies that consecutive step direc- 
tions are correlated [20] • This prediction will be con- 
firmed in a later section, where we analyze stepping data 
from Sowa et al. [l8| for a flagellar motor, and find the 
clear step-step correlations shown in Fig. [5Jb) . In con- 
trast, consecutive cycle completions are statistically in- 
dependent [M Hi. 

The observation that steps arc correlated has impor- 
tant implications for the interpretation of stepping ex- 
periments, and also motivates further theoretical study 
of step ping statistics that goes beyond the assumption 
[H,[l3, 28] of independent steps and dwell times. In this 
paper, we extend the theory of Tsygankov et al. [2(j in 
two ways. First, we consider the distributions of con- 
ditional dwell times instead of mean values. Second, we 
consider a larger class of models, including motor detach- 
ments, substeps, and multiple pathways. 

Our main results are a distribution symmetry for con- 
ditional dwell times, and a simple relation between the 
dissipated free energy per cycle, AG, and certain con- 
ditional stepping probabilities. For a large number of 
models relevant to describe reversible motor proteins, we 
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Here, /9+_|_(i) and p (i) are the probability density func- 
tions for the conditional dwell times between two consec- 
utive forward and backward steps respectively, P±±(t) = 
Jo P±±(t)dt are the corresponding integrated probability 
functions, 7r ++ is the probability that a forward step is 
followed by another forward step, and 7r__ that a back- 
ward step is followed by another backward step. The dif- 
ferent types of dwell times are illustrated in Fig. [TJ and 
Eg. Hal simply means that the conditional dwell times r++ 

and r are random variables with equal distributions. 

However, Eq.[T]docs not say anything about the distribu- 
tions of dwell times between steps of different directions, 
which in general have different distributions. Equation 
[1] holds also at finite average velocity. For example, a 
forward-moving motor will take mostly forward steps, 
but might eventually produce two consecutive backward 

steps. The dwell times between such ++ and step 

pairs have equal distributions, although the probability 
of observing two consecutive backward steps might be 
very small for a motor with strong forward bias. 

In practice, the need to observe a significant number 

of ++ and events to test or apply Eq. [I] can be an 

experimental challenge. For example, ATP driven mo- 
tors like kinesin or myosin V typically have AGatp ~ 

~2hk^T in vivo (29| . This means that events are 

very rare (see Eq. |T8|) . unless an external load is applied 

To extract qualitative information out of an experi- 
mental test of dwell time symmetry, one needs a charac- 
terization of the class of models which satisfy Eq. [T] Wc 
try to formulate a general characterization of this class. 
This means that some of the models within this charac- 
terization will not be realistic descriptions of biological 
systems. On the other hand, a large number of candidate 
models can be excluded if dwell time symmetry is not 
observed in some system. We find two sufficient assump- 
tions for Eq.[T]to hold, which we call strong coupling and 
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Figure 1: Example of the four different conditional dwell times 
in a model with one step per enzymatic cycle (synthetic data) . 
The coordinate x(t) could be position of a linear motor, net 
rotation of a rotary motor, or the net number of a substrate 
or product in an enzymatic reaction, and changes stepwise 
in steps of size d. The staircase line is an idealized running 
average. 



the bottleneck property. Within the discrete state mod- 
eling framework that we use, these assum ptio ns are easy 
to formulate, but can also be disposed of [3fj |. 

The bottleneck property is an assumption about the 
model topology. It means that both forward and back- 
ward steps correspond to transitions to or from a single 
state in the mechanochcmical cycle, which we call the 
bottleneck state. This ensures that the state of the sys- 
tem is uniquely determined after each observed step, in- 
dependent of previous step directions. The bottle-neck 
property alone results in a particularly simple form for 
the step-step correlations, given in Eq. 1151 which can be 
tested experimentally. 

The assumption of strong coupling, defined mathemat- 
ically in Eq. 1201 is related to microscopic reversibility. 
Physically, it means that there is no futile free energy 
dissipation in the mcchanochemical cycle. In addition to 
tight coupling in the usual sense [3fJ, i.e., a one-to-one 
correspondence between fuel consumption and forward 
steps, it includes the assumption that backward steps is 
tightly coupled to synthesis of fuel. As a result, motors 
with strong coupling arc effectively one-dimensional sys- 
tems, in the sense that only one reaction coordinate, e.g., 
position, is needed to describe their operation. 

Strong coupling might not necessarily hold for all mo- 
tors. For example, the forward steps of Myosin V are 
tightly coupled to ATP hydrolysis, while backward steps 
independent of ATP binding was recently reported [8[. 
Another example is kinesin, for which ATP binding dur- 
ing steps in both directions has been reported [l], 0]. 
On the other hand, strong coupling might be possible 
in some rotary motors, at least under certain conditions 

d US HE Ha ■ 

Several recent works [H, [34], [35|, [3(| use Hidden Markov 
models to estimate kinetic parameters directly from ex- 
perimental trajectories, and thus in principle utilize all 
information in the data and bypass the difficulties asso- 
ciated with step detection. However, this approach does 
not replace the need for a theoretical understanding of 
the capabilities and limitations of the underlying stochas- 
tic models. Our results should provide useful guidance 
when applying these techniques. 

Equation [1] corresponds to a similar result in cycle ki- 
netics [H, [2J, [H, [33] • Using the symbol "~" to denote 
cycle time properties, one has 

p + (t)=p-(t), 7T+/*- = J + /J_ = e-^/teT (2) 

where p±(t) are the probability density functions for 
waiting times before forward (+) or backward (— ) cycle 
completions, n± are the relative frequencies of completed 
forward and backward cycles respectively, and J± are 
37 . Equation [5] is based on similar 
24, [25[ . However, the waiting times 



one-way cycle fluxes 
assumptions as Eq.[TJ 
in Eq. [5] do not describe the experimental dwell times 
observed in stepping experiments on motor proteins [20| . 
Distinguishing between cycle completion times and dwell 
times results in significant differences when modeling and 
interpreting stepping data. 
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One example of such a difference is the possibility to 
calculate AG from observations of forward and backward 
steps. In this case, using the cycle completion result in 
Eq. [5] on stepping data can give large systematic errors 
in the estimated AG, which can be avoided if Eq. Q] is 
used instead. 

A related example concerns how to interpret the ra- 
tio of observed forward and backward steps, which has 
been measured for kinesin over a range of forces and 
ATP concentrations [l|, [H, The force dependence 
of this ratio can be described as roughly proportional 
cxp(— cF x d/k B T), where d = 8.2 nm is the step length, 
F x is the applied load, and c is a numerical factor signifi- 
cantly smaller than unity. The dissipated free energy per 
cycle depends on applied load as 



AG(F X ) = E,.d + AG(0). 



(3) 



This means that Eq. [2] predicts c = 1 if steps are assumed 
to correspond to completed cycles. In contrast, Eq. Q] 
gives no general reason to even expect an exponential 
behavior of the ratio of forward and backward steps [2(| , 

and suggests that 7r ++ /ir is a more relevant quantity 

to study. As mentioned above, kinesin might not satisfy 
strong coupling, in which case one must look at more 
complicated models. In any case, making the distinction 
between steps and completed cycles is clearly important 
in order to interpret the experiments correctly. 

The rest of the paper is organized as follows. In the 
next section, we sketch a derivation of our results for the 
models studied in Refs. 0, El]. We also discuss step- 
step correlations. After that, we generalize our results to 
a broader class of models, and discuss detachments and 
multiple pathways in the mechanochemical cycle. We 
also generalize Eq. [T] to the case of several detectable 
substcps. To illustrate our results, we then analyze step- 
ping data for the flagellar motor, before the concluding 
discussion. Detailed derivations are given in the appen- 
dices. 



Sequential Models 

In this section, we derive Eq.[T]and some related results 
in the simplest case, i.e., a sequential model. 

Consider a model with a single cycle consisting of a se- 
quence of N states (see Fig. [5]), where we defined [k]i to 
denote state k in cycle I, associated with position Id. This 
is a basic model in enzyme kinetics [24], (3?J and has been 
used to describe motor proteins like kinesin [23|, [38|, H^] , 
myosin V [22j and Fi-ATPase [40| . By construction, it 
satisfies both strong coupling and the bottleneck prop- 
erty. 

We denote the forward (backward) transition rates 
from state j to adjacent states by Uj (wj) as indicated 
in Fig. [21 The transition rates are positive and periodic, 
i.e., Uj + N = Uj, Wj + N = Wj, and possibly functions of 
external loads and concentrations of various species in 



the surrounding solution. For the purpose of this dis- 
cussion, the rates are assumed to be arbitrary positive 
constants, some of which can be tuned experimentally. 

We assume that transitions between different cycles, 
[TV — 1]; ^ [0];+i, produces observable forward (+) or 
backward (— ) steps. If the model describes an enzyme, 
the observable step could be the release or uptake of a 
product or substrate molecule. 

Counting steps is different from counting cycle comple- 
tions, which can be illustrated with the following thought 
experiment. Consider the sequential model in Fig.r5]with 
N = 2 states per cycle, and a trajectory where the mo- 
tor goes through the states [0]o — ► [l]-i — ► [0]o — * [l]o —* 
[0]i. This completes one forward cycle from [0]o to [0]i, 
but produces three steps: one backward ([0] — > [l]-i), 
followed by two forward ([0]_i — ► [0]o, and [l]o — > [0]i). 
Due to such events, neither the number of steps nor 
the dwell times between steps are accurately described 
in terms of cycle completions. A different treatment is 
needed. 

The time evolution of the system is a random walk 




Figure 2: A simple sequential model of an enzyme or motor 
protein with N states per cycle, [0]i, [2];, . . . , [N — The 
model can be viewed as a biased random walk on the peri- 
odic lattice of states. Numbers in brackets represent different 
states, and indices different cycles. State [k]i denotes the same 
structural state as [k + N]i = but with position or an- 

gular orientation x (for linear and rotary motors respectively) 
shifted by the step length d. Arrows indicate the allowed tran- 
sitions between neighbor states with transition rates Uj, Wj. 
Fat solid arrows denote the major conformational changes 
that produce observable steps like those in Fig. [1] The other 
transitions are indicated by dashed arrows. In a motor pro- 
tein, these transitions could for example be binding or release 
of ligands, or small conformational changes that are hidden 
in the experimental noise. 
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on the periodic one-dimensional lattice of states, where 
the average velocity may be nonzero. To simplify the 
notation, we now use j to denote a state in any cycle, 
with the convention that j and j+N are equivalent states 
in different positions. The probability qj(t) to be in state 
j at time t evolves according to a Master equation [l9j |. 
in this case 

d t qj(t) = Uj-iqj-i(t)+Wj + iqj+i(t)-(uj+Wj)qj(t). (4) 

State j has free energy Gj, and according to detailed 
balance, the free energy difference between two adja- 
cent states are related to the transition rates through 
Uj/wj + \ ~ e~( Gj+1 ~ Gj ^ kBT . From the periodicity of the 
rates we find 

G j+N = Gj + AG, e - AG / k * T = TT — ( g ) 

to m+1 

The direction of the average drift is positive (to the right 
in Fig. \2§ if the dissipated free energy per cycle, AG, is 
negative, and zero if AG = [13, Ijlf . 

Tracking the position produces a scries of forward 
(+) and backward (— ) steps separated by random dwell 
times, as sketched in Fig. [JJ Our aim is to describe the 
statistics of such a trajectory, i.e., the fluctuations of the 
dwell times and the number of forward and backward 
steps, which can then be compared with experimental 
data. 

Following Ref. [2(| , we introduce the pairwise splitting 

probabilities 7r ++ , , n |_, and n , where 7r ++ and 

7r_| is the probability that a forward step is followed 

by a forward step or a backward step, respectively, and 
similar for 7r ++ , tt The splitting probabilities satisfy 

7T ++ + 7T + _ = 7T_ + + 7T = 1. (6) 

Similarly, we introduce random variables t ++ , , t |_, 

and r for the conditional dwell times, where r ++ is 

the dwell time between two consecutive forward steps, 
and so on, as illustrated in Fig. [JJ Once the proper first 
passage problem to describe the dwell times has been for- 
mulated, explicit expressions for the splitting probabili- 
ties and mean conditional dwell times can be computed 
for arbitrary N [20| . using standard methods [l9j . 

The First Passage Problem 

We now describe the first passage problem for dwell 
times in sequential models, which we then use to derive 
Eq. [TJ The dwell time symmetry is due to a one-to-one 

mapping between each event, i.e., two consecutive 

backward steps and the dwell time between them, and a 
corresponding ++ event. This is the key observation for 
the generalization of Eq. [JJ to more complex models in 
later sections. 

We model the splitting probabilities and probability 
distribution functions for the dwell times as a first es- 
cape problem from the interval of states 0, 1, . . . , N — 1. 



One way to approach this problem is the approach with 
absorbing boundaries [l9| . which means solving a re- 
duced Master equation for the states (XI, . . . , N— 1, with 
absorbing boundaries at both ends 0, From the 

assumption that steps are produced by the transition 
[TV — 1]( ?± [0];+i, it follows that immediately after a 
+ step, the system is in the state (±IN), and just after 
a — step, it is in the state N — 1 (±IN). 

We write the reduced Master equation in matrix form, 
dtq = Mq, where the matrix M has elements 

Mij = UjSi, j+ i + Wjdij-i - (uj + Wj)8ij, (7) 

and < i,j < N — 1. The element is the transition 
rate from state j to state i. The probability functions 
are given by the outgoing probability current, 

W±+P± + (t)=U N - 1 [ qN-^dt (8) 
Jo 

K±-P±-(t)=w [ q (t)dt, (9) 
Jo 

where ± indicates an initial condition qj(0), given by the 
direction the of the previous step (see Eq. 159")) . Normal- 
ization of the distributions requires that 

/>oo 

limP±±(t)= / p± ± (t)dt=l. (10) 



Dwell Time Symmetry 

To derive Eq. [JJ we show that the Taylor series of 

7r++P++(£) and tt P__(t) are identical up to a factor 

e AG/k B T^ f rom w hi cn Eq. [TJ follows. The actual calcula- 
tion is given in appendix A, and we now discuss some of 
its consequences. 

First, it is interesting to note that the periodicity of 
the model is not necessary for the dwell time symmetry, 

only to get the simple relation between 7r ++ , tt and 

AG in Eq.[lb] 

Second, the + + / probability distribution has a 

simple closed form for (periodic) sequential models. As 
we show in appendix B, 

N N 

(ii) 

where the Xj are the eigenvalues of the matrix M. This is 
the distribution of a sum A^ of independent exponential 
random variables with mean values AJ~ , , . . . , . 

Explicit expressions for all four conditional dwell time 
distributions for the case N = 3 are given in Ref. [2lj . 
As expected from the calculations in appendix B, there 
is no simple relation between p_| (t) and p \-{t). 

Equation [TTJ is only valid if the eigenvalues Xj are 
distinct. This is the generic situation, since there are 
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no symmetries or other reasons to expect degeneracies. 
However, Eq. Q]is valid also in the degenerate case. 

Finally, it is worth noting that the sequential models 
in this section can also be obtained as a discretization 
of overdamped one-dimensional diffusion [42j in an arbi- 
trary potential U(x) between two points A and B (see 
Fig. [3]). In this case, our results mean that the waiting 
time between last touch at A and first touch at B, and 
the reverse waiting time between last touch at B and first 
touch at A, have equal distributions. This pair of first 
pas sage problems was previously studied by Bier et al. 

Step Directions as a Markov Chain 

The fact that the motor is in a different state imme- 
diately after steps in different directions means that the 
step directions might be correlated. In this section, we 
formulate this observation mathematically, and point out 
some experimentally relevant consequences. 

Let 7r + (fc) and 7r_(fc) = I — 7r + (fc) be the probabili- 
ties that step k in a trajectory is a forward or backward 
step respectively. In the simplest case of no substcps 
or detachments, the definitions of the pairwise splitting 
probabilities give 

/tT+0 + 1) \_ / 7T++ 7T_+ \ ( TT+(k) \ 

{ir-(k + l)J ~ ^7T+_ tt__ J \ir-(k) J ■ [i ^> 

Equation [12] describes the sequence of forward and back- 
ward steps as a Markov chain [H^. I20I] . and the normaliza- 
tion constraints in Eq. [5] leave two independent parame- 
ters in the 2x2 transition matrix. 

The average frequencies tt± of forward and backward 
steps have been measured for several motor proteins [l], 
S B Hi ■ In our model, those frequencies are given by the 
stationary distribution of Eq. [T2] , 




In contrast to the ratio 7r+/-7r_ of the number of for- 
ward and backward cycles in Eq. [21 the ratio ir*^ / n*_ of 
forward and backward steps is in general not equal to 
exp(— AG/fceT). Looking back at Eq. [31 we see no gen- 
eral reason to expect the ratio of forward and backward 
steps to depend exponentially on the applied load. 

Several experiments have divided the dwell times ac- 
cording to the direction of the following step [l], [2, [![. 
Our model gives the forward and backward dwell time 
probability density functions as 

P*+(t) = 7T ++P++ (t) +7T_ + p + _(i), 

p*_(t) = 7r + _p_ + (i)+7r__p__(t). [ > 

Since p^ (t) and p ^(t) are in general not related in 

a simple way, one should not expect equal dwell times 
before forward and backward steps, in contrast to the 
case of forward and backward cycle completion times p3l . 
ii- 

Equation l 1 2l implies that step directions are correlated, 
and gives the step-step correlation function 

C{n) = (s lSl+n ) - ( 8i ) (s l+n ) = C(0) 7 I"I , (15) 

where s m = ±1 indicates the direction of step m, and 

7 = 7T ++ +7T__ — 1 = 1— 7T_| — 7T |_ (16) 

is an eigenvalue of the transition matrix in Eq. [T^] f20| . 
The other eigenvalue is one. 

Note that I7I < 1 since ir±± < 1. In addition, it is rea- 
sonable to expect 7r ++ < 7r |_, since a forward step fol- 
lowing a backward step can be accomplished by a single 
transition, while the whole cycle must be passed before a 
forward step is followed by another forward step. Insert- 
ing 7T++ < 7r |_ in Eq. 1161 and using the normalization 

in Eq. we get 

7 = 7T ++ + 7T - 1 < 7T_+ + 7T -1=0. (17) 

Apparently, we should expect negatively correlated steps, 
— 1 < 7 < 0. Uncorrelated steps (7 = 0) occur in systems 
with N = 1, i.e., simple random walks, as well as some 
special cases of the extended models in the next section. 

An upper bound for the number n of events is 

useful to estimate the regime where our results can be 
applied in practice. A rough estimate is given by 

n__/n tot = tt* tt__ = Trtir ++ e AG / k * T < e AG / k * T . 

(18) 

Hence, the number of steps n tot required to get adequate 
statistics grows at least exponentially with decreasing 
AG. 



Figure 3: The equations of motion for a particle (gray circle) 
diffusing in an arbitrary potential (solid curve) can be dis- 
cretized to a sequential model [12]. In the continuum limit, 
the dwell time symmetry for the discretization translates to 
equal distributions for the two last-touch first-touch waiting 
times between the (arbitrary) positions A and B. 



Extended Models 

Beyond the simplest descriptions in terms of a sequen- 
tial model for the dominating pathway (if there is one), 
more complicated situations with multiple pathways are 
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possible [E BSE SEES SI- A general (coarse-grained) 
description of a motor protein would keep track of both 
the position and the consumption of fuel molecules, which 
requires the model to include an effective chemical coor- 
dinate in addition to the position of the motor (30| . Are 
the dwell time symmetry properties of sequential models 
valid for more general cases as well? As we will see next, 
the answer is yes for a large class of models, in which the 
motion is still effectively one-dimensional. 

With the extended models, we look for general proper- 
ties that imply dwell time symmetry, even if not all mod- 
els with these properties are biologically relevant. We 
find that the strong coupling and bottleneck assumptions 
mentioned in the introduction are sufficient. This means 
that if dwell time symmetry is observed not to hold for a 
particular system, then one can conclude that the strong 
coupling or bottleneck properties are absent. To model 
such a system, one must go beyond the extended models 
presented below, for example along the lines discussed in 
Refs. HE IH, Ho]- The special situation without strong 
coupling but with the bottleneck property is illustrated 
below for an example model (Fig.(7J. 

The extension of our results to detachments and ob- 
servable substeps in later subsections is directly mo- 
tivated by experimental observations of such events 

an. 

The Class of Extended Models 

In this section, we describe a large class of extended 
models which display dwell time symmetry in a little 
more detail, and indicate how the dwell time symmetry 
is derived. We use chemical kinetics, but write a Master 
equation with arbitrary transition rates Wij from state j 
to i, 

d t qi{t) = y~]wijqj(t) -y~]wjjqi(t). (19) 
j¥« 3 

Non-zero "diagonal" rates wa are allowed, for example 
to describe irreversible detachments of motors from their 
tracks [HI, [5l| , or other events that can be filtered out 
exp er iment ally. 

We also need to specify the condition of strong cou- 
pling. One way is to demand that the transition rates 
around any closed loop io, i\, 12, ■ ■ ■ , i m = io of transi- 
tions satisfy 



This property is automatically fulfilled for sequential 
models, and makes it possible to define a free energy Gj 
for each state j, up to an additive constant. From de- 
tailed balance, we have Wij/wji = e~^ Gi ~ G '^ kl>T for the 
free energy difference along the transition i ±=f j. Hence, 
Eq. [201 says that the sum of free energy differences along 
any closed loop is zero, which means that free energy 



dissipation that produces no net motion, e.g., futile ATP 
hydrolysis, is ruled out. (Note that the mechanochemical 
cycle, say, from a state (k) to (k + N), is not a closed 
loop in our description.) We also assume that no tran- 
sition is irreversible, i.e., ^ wji ^ 0, although 
the transition rates may be arbitrarily small. 

To reflect the cyclic operation of motor proteins, the 
transition rates are periodic with period N , and there is 
a well defined free energy per period, 

w 3+Nti+N = Wij, G l+N = Gi + AG. (21) 

An important step in the derivation of Eq. [T]is to con- 
clude that the state of the motor immediately after a step 
is independent of earlier steps. In a sequential model, this 
occurs because there is only one step-producing transi- 
tion per cycle. This condition can be relaxed somewhat. 
It is sufficient to assume the bottleneck property men- 
tioned in the introduction, i.e., that all transitions cor- 
responding to a step are either to or from a single state, 
the bottleneck state, as illustrated in Fig.[4j As the term 
indicates, the system must visit the bottleneck state each 
time it goes through the cycle. 

Whether the bottleneck assumption holds must be de- 
termined for each system separately. For example, devi- 
ations from the step-step correlations predicted in Eq. [15] 
means that the assumption is not valid. 

Although the models considered in this section arc 
more general than the sequential ones, several proper- 
ties of the sequential models are retained. These include 
strong coupling, periodicity, and that steps (roughly) cor- 
respond to one transition in the enzymatic cycle. They 
therefore describe an effectively one-dimensional motion, 
where the motion along the spatial and chemical reac- 
tion coordinates [3(| are tightly coupled to each other for 
motion in both directions. Hence, backward and forward 
motion proceeds in opposite directions along the same 
reaction paths, just as for the sequential models. 

Dwell Time Symmetry and Detachments 

Dwell time symmetry for the extended models with de- 
tachments can be derived using the same methods as for 
sequential models, and further details on that derivation 
are given in appendix C. Just as in the sequential case, 
periodicity is only necessary to establish Eq. Ilbi while 
the dwell time symmetry, Eq. Ilai follows also without 
assuming periodicity. 

Examples of the non-periodic case are the subcycles in 
Fig. which shows a model with two observable sub- 
steps per cycle. Each step in the example satisfies the 
bottleneck property, but both bottleneck states are in 
the same subcycle. Deriving dwell time symmetry for 
the different subcycles proceeds as in appendix C, with 
slightly different initial conditions. 

As mentioned above, detachments of a motor from its 
track are observed in several systems @, SI Bl- W e model 
detachments by introducing death rates [23l.l5l| as extra 
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terms in the diagonal elements of the Master equation 
and M matrix of Eqs. [JJ and SSI 

The presence of detachments affect the long time be- 
havior of the model 5l|, as well as dwell time distri- 
butions and pairwise splitting probabilities. Conditional 
detachment probabilities also have to be added to Eqs. [5] 
and[H 

However, in the derivations of dwell time symmetry in 
appendix C, the transition rates always enter as ratios 
of forward and reverse rates, e.g., Mk m k m+1 /Mk m+1 k m in 
Eq. 1421 so that diagonal elements {k m = k m+ \) always 
cancel. Hence, Eq. Q] remains valid even in the presence 
of detachments. 



Observable Substeps 

For some motors, the full forward and backward steps 
are divided into observable substeps 0, [l4| ■ In this sec- 
tion, we show how the relation between splitting prob- 
abilities in Eq. [1] can be generalized, if strong coupling 
holds and all substeps satisfy the bottleneck property. 

To start with, consider the system in Fig. [5] which pro- 
duces substeps g?i and d,2 — d — d\ during each cycle, and 
assume that each sub-cycle displays dwell time symme- 
try. From now on, we denote a forward step of length 
d\ with +i, and similarly for the other substeps. Keep- 
ing track of only the G?2~steps, we can analyze the system 
as in previous sections, using four splitting probabilities 
7r ± 2 ±2 = n ±± an d dwell times t4 2 ± 2 (£) = T±±(t), which 
satisfy Eq. [JJ 

On the other hand, if we keep track of and discrimi- 
nate between all substeps, we could instead measure eight 
splitting probabilities: 7r +1+2 , ir +1 _ 11 tt_ 1+1) 7r_ I _ 2 , and 
so on. Summing over all paths between two consecutive 



-2-steps, we get 




^0*------"© 

CX v. *© 




© 



Figure 4: Example of a model with the bottleneck prop- 
erty and N = 6 states per cycle. All states inside the gray 
area belong to the same cycle. The bottleneck states are 
0, ±iV, ±2iV .... Black arrows denote transitions coupled to 
a mechanical step of length d, and dashed arrows indicate 
transitions within a cycle, which are undetectable. It is not 
possible to go from one cycle to the next without either leav- 
ing or arriving at a bottleneck state. In general, transitions 
between any states within a cycle are allowed, and bottleneck 
states might have a transition to any state in the neighbor 
cycle. 



+ , 1 + 7T+, 7T_ 



+ 2 + 2 - "+2+! J- T /l+j-i/l-t + i 



7r +2 + i 7A7r+i-i7 r -i+i) fc 7T+i+2 

k=0 



1 - 7T +1 _ 1 7T_ 1+1 



and similarly 



7T_ 2 _ 1 7T_ 1 _ 2 
1 - T+i-iTT-i+i 



(22) 



(23) 



This means that the relation between free energy and 
splitting probabilities for a system with only one step 
per cycle generalizes to 



^+2+1^+1+2 
7r -2-l 71 "-l-2 



_ e -AG/fe B T 



(24) 



Note that we summed over all paths between two +2 
steps without explicit reference to all paths that start 
with a +2 step and end with something else. This means 
that Eq. [24] is valid also for a system with detachments, 
or with several parallel pathways. 

For systems which go through substeps d\ , d 2 ■ ■ • &k m 
each cycle, we can use the same argument to relate the 
splitting probabilities of the full analysis (all substeps 
included) to one where the first substep (d±) is ignored: 



'+K+2 



"+K+1 "+1+2 

1 - TT+i-iTT-i+i ' 

7r -2-i 7r -i-A- 
1 - TT+l-iTT-i+i 



(25) 
(26) 



Iterating this transformation to ignore substeps 
2, 3, . . . , K — 1 as well, we find the following relation for 
a cycle with K visible substeps: 

n +K + i n +i+2 _ _ _ 7r +(g-i)+K _ TT-hg+K _ e -AG/fc B T 
TT-i-K TT-a-i ^-K- {K -i) k'-k-k 

(27) 

In a complicated system with parallel pathways where a 
cycle can be completed using different sequences of sub- 
steps, Eq. [57] holds for every such sequence separately. 





full cycle 



Figure 5: Simple periodic model with two substeps per cycle, 
with step lengths di and tfe. Subcycle 1 has Ni states, and 
subcycle 2 has N — Ni states. Both steps have the bottleneck 
property, although the bottleneck states and Ni — 1 are both 
in subcycle 1. If strong coupling is satisfied, both subcycles 
display dwell time symmetry. 
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An Example Model 

In this section, we consider a small but non-trivial 
model with 2 states per cycle, to present an explicit ex- 
ample that illustrates the results of the previous sections. 

The model has two parallel pathways as sketched in 
Fig. [6l and the steady state velocity and effective disper- 
sion are known exactly (44|. The bottleneck property is 
satisfied, and strong coupling is equivalent to 



UqUiCI = WQWib 



WqWi 



-AG/k-aT 



(28) 



Computing the dwell time distributions and splitting 
probabilities (see appendix D for details), we get 

7T++ = (a(ui + wi) + uqUi)/c , (29) 
7T — = (&(«i + +w w 1 )/c , (30) 

a Aie Al * - A 2 e A2 * 



P++(t) = c 



Ai — A-2 



Ai — A2 



(31) 



e Alt -e A2t b A ie Al * - A 2 e A2t 

P— (t) = C — : r + - r r . (32) 



A1-A2 



Ai-Aa 



We see that p^ (t) and p differ both from each 

other and from the -I — h and distributions, indepen- 
dent of whether Eq. [28] holds or not. The step directions 
are anti-correlated, since 



7 = 7T+H 



— 1 = — WqUi/cq < 0. 



(37) 



The example illustrates the logic of our results. Strong 
coupling and the bottleneck property together are suffi- 
cient conditions for the dwell time symmetry. The bot- 
tleneck property alone is sufficient for the Markov-chain 
description of step directions in Eq. 1121 which remains 
valid even if Eq. [28] is not satisfied. 

A possible reason for why the model in Fig. [6] might 
not satisfy Eq.[28]is illustrated in Fig. [7] This is again an 
N = 2 state model, but with an independent "chemical" 
reaction coordinate (the superscript m on the states), 
e.g, the number of hydrolyzed ATP molecules [3(3]. For 
example, the path [0]]™ ^ ^ [O^j 1 could be a step 
driven ATP hydrolysis, while the reaction [0]" 1 £=; [0]™ j 
could describe an ATP independent step. 

This model can be transformed to the ID model in 
Fig. [6] through the projection = Yl m \j]T> which 



The eigenvalues Ai^ are given in appendix D, and Co = 
A1A2 > 0. For this model, Eq. [25]is obviously equivalent 
to the dwell time symmetry of Eq. [T] For the other two 
pairs of steps, we get 



7T + _ = (b + W Q )(u 1 + W 1 )/c , 

7T-+ = (ui(u + w ) + a(m + W\))/cq, 



(33) 
(34) 



e Alt -e A2t b + w Aie Alt - A 2 e A2 * 
P+-W = c o— : ; 1 : : , (35) 



Ai — A2 



P-+(t) = c-o- 



7T^ — Al — A2 

ab + WQUi Aie Alt - A 2 e A2t 



Ai - A 2 n- + (b + w ) Ai - A 2 



(36) 




u i^^ U ls^' U ]y^^ 
^ y"wo > ywo 

(wH) (wff) (wsD 



w v 



:u 



w v 



;u 



;u 



ui^' U ]/^^ U l^^ 

yy^w^ s ^w® yy^w^ 

CUE) 



W 11 



'WO 



Wv 



;u 



;u 



[ [Q]ri) ^^ nojF) ^==: ( [Q]ftJ 



a 

T 



-d 



Figure 6: Example of a non-sequential model with N — 2 
states per cycle, with two parallel pathways that produce a 
step in the x direction. Solid arrows are step-producing tran- 
sitions, while dashed arrows indicate hidden transitions that 
do not produce an observable step. Shaded boxes are bottle- 
neck states. 



Figure 7: A model with one spatial and one chemical reaction 
coordinate, denoted by indices / and m respectively. Shaded 
boxes are bottleneck states. Projecting out the chemical coor- 
dinate as explained in the text, this model can be brought to 
the same form as in Fig. [fj In general, the resulting effective 
one-dimensional model does not satisfy the strong coupling 
assumption of Eq. [28] 
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leaves the rates unaffected. 



reaction path [0]J" ^ [1]J" ^ [0]™! 1 



However, the open-ended 
[0]™ +1 is trans- 
formed into the closed loop [0]; ±=> [l]j i=> [0]j+i ^ [0];. 
Since strong coupling places no constraints on the rates in 
the original model in Fig. [7j there is no reason to assume 
that Eq. [28] holds for the projected model. In the effec- 
tive ID description, this loop then becomes a slip loop, 
which dissipates free energy without producing net mo- 
tion. Hence, the simple step-step correlations of Eq. [15] 
is retained, but not the dwell time symmetry. 

In systems where it is not possible to measure the 
"chemical position" with single molecule precision, the 
resulting violation of dwell time symmetry is a useful 
test for such projected slip loops, if the bottle-neck prop- 
erty is satisfied. This property can in turn be ruled out, 
for example if the step-step correlations are more com- 
plicated than predicted in Eq. Qj)] 

The most interesting candidates for such tests seem to 
be the rotary motors driven by ion fluxes, e.g., the Fo 
part of ATP synthase, or the bacterial flagellar motor, 
since it is difficult to measure the ion flow with single 
molecule precision. We are not aware of any experiments 
on ATP synthase under conditions that produces both 
forward and backward steps. Data from a flagellar motor 
will be analyzed in the next section. 



Application to the Flagellar Motor 

In this section, we apply our theoretical results to an- 
alyze stepping data from a recent experiment with a chi- 
maeric flagellar motor [la ]. As we will see, the data 
is consistent with the predicted dwell time symmetry 
and step-step correlations within the experimental uncer- 
tainty. We also estimate the free energy per cycle, and 
show that the estimate based on cycle completions can 
have a significant systematic error, compared to Eq. lib! 

The flagellar motor propels many swimming bacteria, 
by driving the rotation of flagellar filaments. Each fila- 
ment is driven at its base by a transmembrane rotary mo- 



ij + 


8 


++ 


+- 


+8 


-+ 




-S 


Trajectory A: 


riij 108 103 


13 


30 


68 


10 


69 


30 


3 


■Kij 0.48 0.46 


0.06 


0.28 


0.63 


0.09 


0.68 


0.29 


0.03 


Trajectory B: 


riij 43 27 


11 


20 


16 


6 


17 


5 


5 


■Kij 0.53 0.33 


0.14 


0.48 


0.38 


0.14 


0.63 


0.19 


0.19 



Table I: Steps (±) and drift (S) events for the two analyzed 
stepping trajectories. A forward step followed by a drift step 
is denoted +8. The rows display the number of events (riij) 
and the corresponding conditional probabilities (jTij). The 
average velocities and average dwell times are zero and 0.04 s 
for trajectory A, and 1.1 Hz and 0.007 s for trajectory B. 



tor, powered by ion flux (Na + in this case [18|]) down an 
electrochemical gradient across the cell membrane. The 
motor is about 45 nm in diameter, and is believed to 
contain 13 torque-generating units [181 l52|. In this ex- 
periment, only one unit was active, and the flagellum is 
expected to switch between 2G distinct orientations per 
turn, corresponding to a step length of 360°/26 w 14° 
fl^ . We analyzed one trajectory with close to zero net 
velocity, and one with finite velocity. 

Aided by Chung-Kennedy filtering [53| to enhance 
the steps, and the step-finding algorithm described in 
Ref. [la] , the raw data was converted to stair-case step- 
ping data, as shown in Fig. [SJa). Upon close inspection, 
the trajectory seems to be divided into intervals where 
the step lengths are consistent with each other, but out 
of phase with steps outside. This apparent drift mig ht be 



54[ that 



due to dynamical exchange of the stator units 
anchor the motor to the cell. 

We identified intervals of consistent stepping by in- 
spection (see Fig.[5Ja)), and excluded the drift from our 
analysis by treating inconsistent steps as detachments, 
i.e., each interval of consistent stepping was treated as 
an independent run. The result is summarized in Tab. [I] 
and we now proceed to compare step-step correlations 
and dwell time distributions from trajectory A with the 
theoretical predictions. Trajectory B contained too few 
steps to make such comparisons meaningful. 

The step-step autocorrelation function for trajectory A 
is shown in Fig.^b). The prediction C(n)/C(0) = 7™ of 
Eq. is valid for trajectories without detachments. In 
order to compare theory and experiments, the splitting 
probabilities in Tab. [J must therefore be renormalized by 
a factor (1 — n±s)~ l to account for the fact that inter- 
vals of consistent stepping (by definition) contain no drift 
events. This gives the theoretical prediction 



7 



'++ 



7T+.5 



1 - TT-a 



- 1 = -0.43 ±0.06. 



(38) 



A least squares fit to the experimental correlations gives 
7 = —0.40, in agreement with Eq. [35] but inconsistent 
with uncorrelated step directions. 

Having confirmed that the step-step correlations are 
consistent with Eq. [T5] we go on to test the dwell time 
symmetry. The empirical distribution functions P±±{t) 
are plotted in Fig. [Hfc) for trajectory A. Error bars are 
standard deviations from bootstrap estimates [55| . Us- 
ing the Kolmogorov-Smirnov test [56j, we conclude, with 

95% confidence, that t |_ have different distribution than 

t++ and t , and that r_| have different distribution 

than t However, the test could not detect significant 

differences between the other three pairs (+ + / , 

H / , and H / — h). This is again consistent 

with the theoretical predictions, and also shows that the 
statistics is good enough to detect differences between 
distributions that are not equal. It is also interesting to 

note from Fig.[5](c) that r ++ and r__ arc longer than r H 

and t |_ on average. This is reasonable, since H — / — h 

cvents in principle only requires a single transition, while 
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the system must go through a complete cycle to complete 

a ++ or event. 

The ionic motive force was not measured indepen- 
dently [3, so it is not possible to compare the free- 
energy estimate of Eq. [1] with an experimental value. 
However, we can compare our result with that based on 
Eq. [2j which was used in Ref. [HI ■ The results are sum- 
marized in Tab. [TTJ As expected, the free energy per 
cycle is close to zero in trajectory A, which has almost 
no net velocity, while trajectory B clearly has a finite free 
energy per cycle to drive the rotation. Also note the dif- 
ferences in the estimates based on Eq.[T]and Eq.[^l which 




n [steps] 



(c) -, 



sr 0.5 




Figure 8: Analysis of stepping data from trajectory A. (a) 
Stepping data and result of step detection (staircase line). 
Grid lines indicate the theoretical step length d = 13.8° . Step- 
ping is consistent with the theoretical step length in intervals 
interrupted by drift events (arrows), and adjacent intervals 
are out of phase with each other, (b) Step-step autocorrela- 
tion function as defined in Eq. 1151 The measured correlation 
function (squares) is consistent with the theoretical prediction 
(circles) using conditional splitting probabilities from Tab. [T] 
Symbols are the same size as the estimated standard devia- 
tions, (c) Distribution functions P±±(t) for the four different 
dwell times. 



show that the cycle-completion estimate can give rise to 
significant systematic errors in estimated free energies if 
applied to stepping data. 



Conclusion and Outlook 

The statistical properties of steps in reversible molecu- 
lar motors are not the same as that of cycle completions 
in enzyme kinetics. To interpret stepping trajectories 
correctly, this difference must be taken into account. 

In this paper, we have extended the theory for dwell 
times in sequential models [2(], [2l| , and derived a symme- 
try relation for the conditional dwell time distributions 
for a large class of experimentally relevant models, in- 
cluding parallel pathways, visible substeps, and detach- 
ments. In contrast to the statistics of cycle completions 
[23l . [24], [HI, the dwell time distributions and splitting 
probabilities of steps depend on the step directions of 
both steps before and after the dwell period. 

The dwell time symmetry is a consequence of strong 
coupling, i.e., tight coupling for both forward and back- 
ward steps, and a bottleneck property of the underlying 
kinetic scheme. This means that the motion of the sys- 
tem is essentially a one-dimensional random walk in the 
(two-dimensional) space of spatial and chemical reaction 
coordinates. This makes it possible to use our results to 
infer information about the coupling and efficiency of a 
system from kinetic data, without detailed assumptions 
about the underlying kinetic scheme. In this respect, 
it is similar to the relation r > 1/N between the ran- 
domness parameter r and the number of steps N in the 
mechano chemical cycle of molecular motors [28l . [57], [58| . 

An alternative to the discrete formalism used here, is 
to include continuous spatial degrees of freedom, as is 
common for example when modeling ratchet-type motor 
mechanisms [3(1 O, H(J ■ Since such models can be dis- 
cretized [42|], it is in principle possible to inquire about 
strong coupling and the bottleneck property in continu- 
ous models as well. An example of this is diffusion in one 
dimension (see Fig. which strictly satisfies strong cou- 
pling and the bottleneck property. An interesting prob- 
lem for further research is to quantify how well the dwell 
time symmetry is preserved in system with small devia- 
tions from these assumptions. 

We analyzed stepping data from single motor subunits 
in a flagellar motor [18[ . The data seems to be consistent 



AG; 



steps 



AG, 



cycles 



Trajectory A 
Trajectory B 



0.06±0.23 -0.05±0.10 
-0.94±0.28 -0.47±0.13 



Table II: Free energy estimates for trajectories A and B, in 
units of ksT. The estimate AG s t cp s is based on our main 
result, Eq. [1] while AG cyc i cs comes from Eq. [2] i.e., treating 
steps as cycle completions. 



11 



with the predicted dwell time symmetry and step-step 
correlations, although only one trajectory with almost 
zero velocity had enough steps to make such a comparison 
meaningful. 

The form of the step-step correlations is consistent 
with an underlying kinetic scheme that satisfies the bot- 
tleneck property So, does the dwell time symmetry 
in Fig. [S] indicate that the flagellar motor subunits are 
strongly coupled to the driving ion flow? Another possi- 
bility is that the steps are equilibrium fluctuations, as in- 
dicated by the low AG estimated from the stepping prob- 
abilities. In equilibrium, there is no free energy change 
associated with ion transport, which means that dwell 
time symmetry could be obtained also for a system with 
loose coupling. To say anything about the coupling in the 
flagellar motor, the dwell time symmetry must be tested 
in a regime with finite (and constant) velocity and ionic 
driving force. This is in principle a question of observing 
more steps in such conditions, but a systematic way to 
identify and separate steps from drift events would also 
be useful. 

The step-step correlations predicted in Eq. [TO] might 
be present also in systems that violates the dwell time 
symmetry. Deviations from this form would say some- 
thing about the topology of the kinetic pathways in such 
systems, and it would probably be useful to apply our 
analysis also to kinesin and myosin V. 

An accurate theoretical model is often crucial in or- 
der to correctly interpret experiments on such complex 
systems as motor proteins. Many previous theoretical 
works on steps and dwell times in molecular motors 
[22l . I23I , l26l . [27j derive or assume descriptions where con- 
secutive steps or cycles are statistically independent of 
each other. As we have demonstrated, this does not ap- 
ply to the flagellar motor of Ref. [l|| . We have also pre- 
sented a large class of simple models where easily acces- 
sible quantities like the dwell times and step directions 
are correlated. 

We expect our results to be of practical use in both 
data analysis and design of experiments. In particular, 
further experimental efforts are motivated in order to de- 
tect correlations and collect significant statistics for both 
forward and backward steps. 



APPENDIX A 

We now give a detailed derivation of Eq.[T]for a sequen- 
tial model. The initial condition is the state immediately 
after a step, and can be written 



qj(0) = = Sjfi, after a forward step, and 
qj(0) = Qj = Sj^—i, after a backward step. 



Eq. [3S This gives 

n ++ d t n P ++ (0) = u N ^i(Al n ~ 1 q < - + ^) N _ ± 

= u N „ 1 ^2 M N - 1>kn _ 2 ■ ■ ■ M k2tkl M klfi , (40) 

where the summation over ki, fe, • • ■ goes from to N— 1. 
The same calculation for P__(4) yields 

n^d t n P—(0) = w {M n ' 1 q^) Q = 

= M okl M klM ■ • • M fc „_ 2 ,jv-iWo • (41) 

{kj} 



The products of matrix elements in Eqs. (|40II41[) corre- 
spond to n-step paths between states and N — 1, plus 
the extra escape step. For sequential models, the short- 
est such path is N steps, so the first N — 1 derivatives 
are zero. For n > N, we note that there is a one-to-one 
correspondence between the nonzero terms in Eq. 1401 and 
Eq. [HJ For each term including a path from to N — 1 in 
Eq. 1401 there is a corresponding term for the reverse path 
from TV — 1 to in Eq. [IT] The ratio of two corresponding 
non-zero terms is 



R({kj}) 



M kl<0 M k2M 



M, 



JV-l,fc„_ 2 UN-1 



M 0>kl M klM M kn _ 2iN -i w 



(42) 



Going back to Eq. [7] for the elements of M, and us- 
ing wq = wn by periodicity, we see that R({kj}) = 
e - AG / k BT f or a vj p a j rs f corresponding terms. More- 
over, P_|__|_(0) = P__(0) = 0, as there is no transition 
which produces two steps at once. Therefore, 

n ++ d t n P ++ (0) = e - AG / fcBT 7r„a t "P_(0) (43) 

for all n > 0. The underlying reason is that the sum of 
free energy changes is the same along all possible paths 
going forward one cycle from 0. 

Since P++(i) and P__(t) are integrals of the solutions 
of the reduced Master equation, which is a finite sys- 
tem of ordinary differential equations with constant co- 
efficients, they are smooth functions which have Taylor 
series. Hence, Eq. 1431 together with the normalization in 
Eq. [TO] implies Eq. [TJ A similar correspondence between 
paths was used in Ref. [25| to derive Eq. [5J The period- 
icity is necessary to get the simple relationship between 
7r ++ , 7r__ and AG in Eq. Ilbl but not to establish the 
dwell time symmetry in Eq. Hal Without the periodicity, 
Eq. [Tb] is replaced by 



= R({kj}) 



= MjV - 1 c -(G W - 1 ~G )/fc B T^ 
Wo 



(44) 



(39) 



APPENDIX B 



To derive Eq. [1] we compute 7r ++ 9 t "P ++ (0) and 
7r__(9 t "P__(0) using Eq. [5] and the initial conditions in 



To derive explicit expressions for p ++ (t) and p (t) 

for sequential models, we first note that the occupation 
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probabilities qj(t) in Eq.[8]are the solution of a system of 
linear ordinary differential equations with constant coef- 
ficients, namely the elements of M in Eq. [7] If the eigen- 
values {Xj} of M are non-degenerate, the solutions have 

the form qj(t) = J2iLi a [^ eXit - Note that the eigenval- 
ues have negative real parts, to guarantee that the system 
eventually leaves the interval. We make the ansatz 



N 



P+ + (t)=P__(i) = l + ^> 



A,/ 



(45) 



where the term 1 ensures proper normalization. We need 
./V equations to determine the coefficients otj . These are 



d t k P ++ (0)=d t k P-(Q) = Q 



(46) 



for < k < N — 1. The case k = was argued just 
before Eq. g31 and 1 < k < N — 1 follows from the 
argument between Eqs.l4"Tl and |4"2"1 The resulting system 
of equations for otj is of the Vandermonde type, 








/ 1 

Ai 
A? 


i 

A 2 • 
Ai • 


' 1 ^ 

Atv 
X 2 


( ai ^ 
a 2 


Ur 1 


\N-1 
A 2 


■ A N 1 


\a N J 



V o J 



Solving with Cramer's rule and setting 



iY 



p++(t) = p—{t) = dtP—{t) = ^2x k 



"fee 



(47) 



(48) 



k=l 



we arrive at Eq.QT] The distributions of t_| and r_| can 

be computed with the same ansatz, but the derivatives 

<9 t fe P_l (0) and d t k P^ (0) are non-zero for k > 1, so the 

results are more complicated and there is in general no 
simple relation between p_| (£) and p \-(t). 



APPENDIX C 

In this appendix, we derive Eq. [1] for a periodic ex- 
tended model like the one in Fig. 2J generalizing the first 
passage problem as we go along. The derivation proceeds 
much as for the sequential model, by solving a reduced 
Master equation for the states 0, 1, . . . , N — 1, 



d t q=Mq, M v 



E 



(49) 



with absorbing boundaries at —1 and N, and allowing 
wu > to model detachments. Note that the sum in 
Eq . |4"91 have at most 2N — 1 terms, due to the bottleneck 
property. After a forward step, the system is in state 0, 
while it could in principle be anywhere in the cycle after 



a backward step. The initial conditions, describing the 
distribution of states just after a ± step, are therefore 

N-i 

qM = q { + ] = iM = lY ] = ^> Z ~ = E WkN - 

k=0 

(50) 

Since there are several escape transitions for each step, 
Eq. [5] for the dwell time probability functions is general- 
ized to 



N-l ft 

■n±+P±+{t) = V w Nk / q k (t)dt 
k=o Ja 

,t -i 
Ti±-P±-{t) = z'_ I q (t)dt, z'_ = w k0 . 

JO ,,_ ,*r 



(51) 



k=-N 



We now proceed to compute derivatives of P++(0) and 

P (0), and note that by Eqs. (|20H21|) , the free energy 

difference Gn — Go = AG is independent of path: 



JV-l 



n ++ d t ' 1 P ++ (0) = ^^{M^qM)^ 

■M klko 6 ko ,o, (52) 



k n -i=0 



E w Nk n - 1 M kn _ lkn 
{k 3 } 



7r__a t "P-(0) = z_(M"- 1 ( f(-)) 

= V M akl M klk2 . . . M kn _ 2kn _, Wkn - lN 

^ 7. 



{ kj } 



= f^ e (G«-Go)/fc B T a „p ++(0) _ (53) 



From the periodicity of Wij expressed in Eq. [2T1 we get 
z'_ = which proves Eq. [33] and thereby Eq. [1] The 
derivation for the case when backward steps (instead of 
forward steps) end in the bottleneck state is analogous. 

Note that our definition of the bottleneck property 
means that there are no transitions that produce two 
steps at the same time. This rules out discontinuities at 
t = in the distribution functions P±±(t), and justifies 
our use of Taylor expansions. 



APPENDIX D 

In this appendix, we solve the example model in Fig.[S] 
in some detail. In addition to deriving Eqs. (|29ti37[) . this 
also illustrates a method that can easily be generalized 
to larger systems and implemented on computer algebra 
systems. For N > 5 states per cycle, the eigenvalues 
Xj must be calculated numerically. With this exception, 
the method gives explicit expressions for the dwell time 
distributions in the real time domain, that can be directly 
compared to experimental dwell time histograms. 
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The first escape problem for the model in Fig. [6] is 
governed by a reduced Master equation as described in 
Eq. I49| with the matrix 



M 



-(a + b + u + w ) wi 

u —(ui + wi) 



(54) 



The eigenvalues of M are 



Al,2 — 



—c\ ± \J (wq + u - u>i - ui + a + b) 2 + 4:U wi 



and satisfy the relations 



(55) 



c = AiA 2 = (wq + a + b){ui + w\) + uom, (56) 
ci = — (Ai + A2) = a + b + uo + wq + m + w\. (57) 

The initial conditions after forward and backward steps 
are given by Eq. 1501 and the conditional dwell time dis- 
tributions are given by Eq. 1511 

*j +) (0) gj- ) (0) = Mo i + ^ > (58) 



b + Wq 



7T±_P(t)±- = (w Q +b) / q (t')dt', 
Jo 

K±+P(t)±+ = / aq (t') + u\qi(t')d£ '. 
Jo 



(59) 
(60) 



Using the ansatz in Eq. [15] together with Eqs. (|58H60p 
to compute dtTt±-P±-{t) and dtTt±+P±+(t), we get the 
following systems of linear equations: 



Finally, we need expressions for the splitting probabili- 
ties, which we derive using the adjoint equations [l!^ | . Let 
7Tj±, with j = 0, 1, be the probability that a system start- 
ing in state j will next produce a ± step. Starting with 
7Tq-, we note that from state 0, a backward step can be 
accomplished either by the next transition being a back- 
ward step, which has probability (b+Wo) / (a+b+Uo+Wo) , 
or by jumping to state 1 and from there eventually get a 
backward step. The last possibility has total probability 
tti-Uq/ (a + b + uq + wq), so that 



7T - 



b + w 
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a + b + uq + wq a + b + u + w 



7T1- 



(62) 



Next, we multiply with ^ fc wok = a + b + wq + uq and 
rearrange the terms. Applying the same procedure for 
the other tt,-±, we finally get 



M 1 



( TTO- 


7T0+ ^ 




\ 7Tl_ 


7T1+ J 





(63) 



b + wo a \ 
ui J 

To get the pairwise splitting probabilities in Eqs. (|29H37p . 
we first solve for the ttj± , and then weight them according 
to the initial conditions of Eq. [58j 



7T ±+ = 2_^ n 



3 + 



(±) 



(0). (64) 



Acknowledgments 




aV a, a. 



a, a.n a 9 a g 



-1 -1 -1 -1 

a b-\-WQ ba-\-WQUi b 

7T ++ 7T + _ (b-\-Wo)7T 1. 7T__ 



(61) 



Solving this and using p±±(t) = dtP±±(t), we eventually 
arrive at the probability density functions given in Eqs. 
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